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ABSTRACT 

We present Keck/NIRC2 K t band high-contrast ooronagraphic imaging of the luminous debris diak 
around the nearby, young A star HD 32297 resolved at a projected separation of r = 0.3-2. 5" (*» 35- 
250 AU). The diak is highly warped to the north aod exhibits a complex, “wavy" surface brightness 
profile interior torn 110 AU, where the peaks /plateaus in the profile* are shifted between the NE 
and SW dink lobes. The SW side of the disk is 50-100% brighter at r « 35-50 AU, and the location of 
its peak brightness roughly coincides with the disk's mm emission peak. Spectral energy distribution 
modeling suggests that HD 32297 has at least two dust populations that may originate from two 
separate belts likely at different locations, possibly at distances coinciding with tbe surface brightness 
peaks. A disk model for a single dust belt including a phase function with two components and a 5 10 
AU pericecter offset explains the disk’s warped structure and reproduces some of the surface brightness 
profile's shape (s.g. the overall “wavy" profile, the SB peak/ plateau Bhifla) but more poorly reproduces 
the disk’s brightness asymmetry and the profile at wider separations (r > 110 AU). Although there 
may be alternate explanations, agreement !>etween the SW disk brightness peak and disk’s peak mm 
emission is consistent with an over density of very small, sub-blowout- siaed dust and large, 0.1-1 mm- 
aized grains at % 45 AU tracing the same parent population of plonctesimals. New near-IR and 
sub in m observations may be able to clarify whether even more complex grain scattering properties or 
dynamical sculpting by an unseen planet are required to explain HD 32297’s disk structure. 

Subject headings: stars: early- type, planetary systems, stars: individual HD 32297 


1. INTRODUCTION 

Debris disks are signposts of planets and planet forma- 
tion (e.g. Wyatt 2008; Kenyon and Bromley 2008). Sup- 
port.ng this picture, the two stars with independently 
confirmed, directly Imaged planetary systems HR 8799 
and 0 Pictoris (Marais et al. 2008, 2010; Currie et al 
2011a; Lagrange et al. 2010) are surrounded by luminous 
debrs disks (Smith and Tbrrile 1984; Rbee et &1. 2007; 
Su et al 2009). Similarly, Fbmalhaut has a candidate 
planetary companion located just interior to the star's 
bright debris ring (Kalas et al. 2008). 

In the absence of a directly imaged planet, re- 
solved imaging of debris disks may provide indirect ev- 
idence for a massive planet’s existence, may holp con- 
strain the unseen planet's properties, and thus can 
help identify promising targets for future direct imag- 
ing (e.g. Wyatt et al. 1999). For example, tbe in- 
clined or “warped” component of 0 Pictoris ’s debris 
disk 'Heap et al. 2000; Golimowskl et al. 2006) is likely 
due to the directly imaged planet (Augercau et al. 2001; 
Dawson et al. 2011) and also provides an estimate for 
the planet's tnasB independent of planet cooling mod- 
els (Lagrange et al 2010). Dynamical sculpting by a 
planet/planets may explain the sharp inner edge and 
percenter offset of Fomalhaut’a debris ring (Kalas et al 
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2005; Quillen et al 2006; Kalas et al 2008). Other de- 
bris dnk structures may be due to non-planet processes, 
in particular interactions with tbe interstellar medium 
or perturbations from a nearby star, as has been pro- 
posed to explain images of disks around HD 15115 and 
HD 61005 (e.g. Kalaa et al. 2007; Hines et al. 2007). 

The nearby (d = 112t}J pc; van Lccuwen 2007) A5 
star HD 32297 is another example of a young star sur- 
rounded by a luminous, spatially- resolved, debriB disk. 
At 30 Myr old (Kalas et al. 2005), it is roughly co- 
eval with HR 8799 and may probe debris disk evolution 
at a stage just after they are most colhsionally active 
(Kenyon and Bromley 2008; Currie et al. 2008, 2009). 
Like 0 Pic and HR 8799, HD 32297 has s large infrared 
(IR) excess emission due to circumstellar dust first iden- 
tified from IRAS data. Schneider et al. (2005) thus se- 
lected HD 32297 for Hubble Space Telucope{HST) NIC- 
MOS (Fl 10W) coronagraphic imaging and resolved the 
disk out to an angular distance (from the star) of ~ 3 3" 
(~ 400 AU). HD 32297 was subsequently resolved in the 
optical (Kalas 2005), near-IR (1.6- 2 2 pm, Debes et al. 
2009; Mamet et al. 2009), thermal infrared (10 20 pm, 
Moorhen et al. 2007; Fitzgerald et al. 2007), and mil 
li meter (1.3ram Mancea et al. 2008L 

Previous work has claimed that HD 32297’s disk struc- 
ture is shaped by planet sculpting as well as non-planet 
processes. Debes et al. (2009) identified an asymmetry 
in tbe disk scattering efficiency between the northeast 
and southwest sides (see also Kalas 2005). They ar 
gued that ISM sculpting explains this feature much like 
it explains some properties of the HD 15115 and HD 
61005 disks Schneider et al. (2005) identified a bright- 
ness asymmetry between the two disk sides, a feature 
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consistent with sculpting by a miaivi planet (m» also 
Mi'iM e*. aL 2008). 

Tne two mechanisms. ISM sculpting and planets, are 
not mutually exclusive. New images of the HD 15116 
and HD 61005 disks reveal cleared inner regions and/or 
per center offsets, both of which are plausibly due to a 
platctary companion (Kodigas et al 2012; Buenxk et al. 
2010). Additionally, multiple debns belta, scaled ana 
logics to the solar system’s asteroid belt and Kuiper beit, 
are also likely planet signposts and may reside around 
HD 61005, HD 15115, and HD 32297 (FiUgerald et al. 
2010; Martess et al. 2008; Rodigas et al. 2012). 

Tb determine which mechanisms are responsible for 
shaping HD 32297’s debris disk structure, we need new, 
high signal- to- noise images with which to derive precise 
disk properties. Although Schneider et al. (2005) kfon 
tify a disk brightness asymmetry consistent with planet 
sculpt ins. they caution that the disk brightness measure- 
ments dose to the coronagraphic spot (r 0 3-0. 4") 

which provide the basis for this asymmetry are highly 
uncertain PSF subtraction errors due to the completely 
opaque NICMOS coronagraphic spot may limit our abil- 
ity to conclusively identify disk structure at these small, 
spedde-dominaLed separations. Moreover, if the ssytn- 
identified s planet-induced density structure, tt 
align with the mm emission peak (Maneas et al. 
2006). However, it is not clear whether these asymme- 
tries are aligned and thus whether they identify small 
and large grains originating from the same parent popu- 
lation The Paloma x/K, Image from Mawet et al. (2009) 
lias limited spatial resolution. While they did recover 
Schneider et al ’s brightness asymmetry, higher spatial 
reso.ution observations could confirm and help clarify the 
physical origin of this and other features. Fbr example, 
new data could identify breaks in the disk bri gh t pees pro- 
file that may reveal evidence for the multiple debris belts 
inferred from modeling unresolved IK data. 

Tb further clarify the nature of the HD 32297 debrw 
disk, we present new coronagraphic imaging obtained at 
K, t~ 2 pm) with the Keck telescope on Maims Kra. 
resolving the disk a: angular separations of 0 3-2 5". \2 
describes our observations and extraction of tbs disk iro 
ages In (3, we investigate basic disk properties (position 
angle, full-width half-maximum, and surface brightness) 
as a function of angular separation from the star We 
then combine imaging with new, unresolved broadband 
photometry from the Spitzer Space Telescope and tbs 
WISE satellite to constrain the disk structure and idnn 
tify the loeation(s) of the disk emission ($4) Finally, 
we compare our analyses to those from previous work on 
HD 32297 ((5) and investigate the physical mechanisms 
responsible lor sculpting the disk emission ($6). 

2. NIRC3 DATA 
2.1. Observations 

We imaged HD 32297 on UT November 20, 2011 with 
the N1RC2 camera mounted on Keck II using tbs K, fil- 
ter (A. • 216 pm) in the narrow camera mods (9 952 
mas /pixel Yelda et al. 2010) with correlated double sam- 
pling The Keck AO system delivered diffraction limited 
images with a FWHM of ~ 4 9 pixels (~ 49 mas). Tb en- 
hance our ability to extract the HD 32297 disk emission 
from the bright stellar halo, we centered the star behind 


the 0.6" diameter, partially- traasmasive coronagraphic 
spot and used the ■ large hex" pupil plane mask 

Our HD 32297 data consist of coadded 30s exposures 
with a cumulative integration time of 20 minutes and 
wore taken through transit (HA » (-0 22,0.36]) in “ver- 
tical angle” or angular differential imaging mode (ADI, 
Marois et al. 2006). Over the course of those exposures, 
the parallactic angle changed by 36 degrees. While light 
cirrus caused fluctuations in the source transmission (as 
measured from the PSF core twhind the coronagraphic 
mask) on the order of ~ 10%, observing conditions were 
otherwise stable. 

Tb flux calibrate the data, we measure the flux of the 
PSF core of HD 32297 as viewe d through the partially 
transmissive coronagraphic mask and corrected for the 
extinction through the mask Tb determine the extinc- 
tion, we use Keck/NIRC2 K, observations of HD 15115 
taken in August 2011 in photometric conditions with the 
same coronagraphic spot stae that were flux-calibrated 
with a standard star HD 3029 Rum the August data, we 
measure the extinction through the 0 6" spot to be 6.91 ± 
0. 15 mags, sagbtly leas than but comparable to extinction 
estimates in A', for 1 2" spot diameters (Metchev at al 
2009). Uncertainties in the aperture correction Inside 
the coronagraphic spot due to scattered light from the 
spot edge limit the precision of our estimate We did 
observe a photometric standard (p486r) 90 minutes after 
our HD 32297 imaging sequence Although conditions 
became highly variable during the standard star obser- 
vations due to patchy cloud coverage, our abeolule flux 
calibration Implied from the least extincted standard star 
frames agrees with that derived from the coronagraph 
transmission to within 0.3 magnitudes. 

2.2. Image Processing 

Basic image processing and high contrast imaging 
techniques follow methods outlined In Currie et al. 
(2010, 2011a,b). We employ standard dark subtraction 
and flatfiekkag c or rec t ions. Identify and interpolate over 
hot/cold pixels, and apply the distortion correction de- 
termined from Yelda et al. (2010). Fbr image registra 
lion, we exploit the fact that the PSF core Is visible 
through the coronagraphic mask and is unaaiurated. We 
register each image to subpixel [o (mr m 0.1 pixels) ac- 
curacy by determining the centroid position of the first 
image in the sequence and then determining the rela- 
tive offset between the reference image and subsequent 
images by solving for the peak in the cross-correlation 
function fbr each image pair. Finally, we subtract off a 
radial profile to remove the smooth seeing halo and make 
a second pass through the images to identify remaining 
bad pixels. 

Tb extract a detection of the HD 32297 disk, we adopt 
the Locally Optimized Combination oj Images approach 
(LOCI Lafreniere et al. 2007), using an updated ver- 
sion of the LOCl-baaed code employed in Currie et al 
(20Ub), which will be detailed later (Currie et al. 2012, 
T Cume, 2012 m prep ) Following Thalmann et al 
(2011), we reduce the data using “conservative" LOCI 
settings more appropriate for extended sources (i.e. 
(inks, not planets), where we oonsidar rotation gaps 6 of 
15-5 x FWHM, optimization areas (N*) of 1 000-3000 x 
the FWHM area, optimization geometries g of 0.5-2, 
and subtraction annuli (dr) of 5-10 pixels wide (see 
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Lafrmbre et al 2007, for more detub), lb determine 
the eifnal- to- notec per resolution dement of our dbk de- 
tect too, we convolve the image with a boom use equal to 
** 1 FWHM and compare pixel counts to the standard 
deration of counts w-.thin a 1 FWHM- wide annulus 10 
Fin illy, we correct lor photomelnc/astrometrtc biaaea in- 
bernnt in LOCI -based processing by unputlnj{ fake disk* 
intc each registered image and comparing the predicted 
and measured disk properties (e.g. surface brightness, 
full- width half maximum, position angle). Our method 
follows that first developed by Rodigaa et al. (2012) and 
b described In full in the Appendix 

2.3. HD nt97 Disk Image 

F.gur« 1 abuwrs the image and signal- to^ noise map for 
our ‘conservative* LOCI reduction «mming a rotation 
gap criterion of 6 £ 2.5, N A ~ 3000. g - 1. and dr m 10, 
which balances our ability to detect the disk at small sep- 
arations by attenuating speckles (favoring smaller 4, S' A ) 
but without oversubtracting the dbk (favoring larger \ 
N a ) As evidenced by the signal- to- none map. we detect 
the disk at SNR > 3 from r =» 2 5" ahthr way to the 
edge of the coronagraphic spot at 0 3" The peak wgnal- 
to-nobe per pixel b — 19. The disk emission along the 
mid plane b more than 10o significant at 0 85-1.4", while 
some regions of dbk emission at 0.& 0 5" on the SW side 
are itill more than 5-<r significant. Visually inspecting 
tha image and signal -to- noise map reveals significant dbk 
alructure Most notably, the disk emission traces a dis- 
tinct “bow” shape, where the disk position angle clearly 
changes with radial separation. 

Furthermore, the disk exhibits a significant brightness 
asymmetry at small separations (r < 0.7"). Figure 2 
redisplays our dbk Image with a different color stretch 
to better illustrate the brightness differences between the 
two rides of the dbk. Interior tors 035" (identified 
with a arcb), the SW aide s significantly brighter than 
both the NK ride as well as any pixel value for residual 
specklea located at different azimuthal separation., the 
NE ride has a peak brightness only slightly larger than 
tbs brightest spedtb Exterior to this separation, there 
are so residual speckles aa bright as either side of the 
disk, and tbs SW b still clearly brighter at least out to 
r os 0 8" (yellow/red region on the SW ride). 

Our Keck K, image agrees well with the previous 
best K, results, which were obtained with the extreme- 
AO Well- Corrected Sub- Aperture on Palotnar presented 
by Mnwet et al. (2009) using classical (not ADI) imag- 
ing Convolved to the PaJom&r/WCS beamsbc (Figure 
2, lx>’.tom panel), our image strongly resembles that of 
Mawet si a).'*. The dbk appears highly asymmetric with 
the SW side being brighter extending all the way to the 
oorocagraphlc spot (0.3" in our unage compared to 0 4" 
in theirs) As with the Mawet et al. (2009) data, the NE 
ride appears fainter and truncated. Moreover, the bright- 
est pevtioo of the dbk on the SW side roughly overlaps 
with the mm continuum peak (Maneas et al 2008). 

3. ANALYSIS 

“The btw at Oils pap«r u taafxcg tad dharaclavtsuw 'J» HD 
1229T VhrW dak While «r> do not pnseat daU iwductloM with 
SMlhocs mare optanfcsd lor point source detection, we pie* to do 
so la e future work 


Here, we investigate the HD 32297 disk geometry 
and surface brightness profile, following methods smu- 
j* to those described in Rodigaa et al. (2012). We per- 
form analysis on the image shown in Figure 1. The 
disk properties we report are corrected for photo met- 
ric/aatro metric biases inherent in LOCI processing (see 
Appendix). 

3.1. Duk FWHM 

Tb better assess the HD 32297 dbk morphology, we 
measure the disk FWHM perpendicular to the dak's ma- 
jor axis aa a function of stellocentric distance. First, we 
identify the brightest pixels at each radial separation for 
the NE and SW lobes, respectively. Next, we place a 6 
pixel by 21 pixel box centered on the brightest pixel and 
sum up the counts /pixel along each row of the box, pro- 
ducing a ID array of 21 values. Finally, we fit a Gaussian 
to this array, which yields the dbk mid plane location and 
the dbk FWHM at that location 

Figure 3 displays the dbk FWHM aa a function of stel- 
locentric <Lstance for the NE (purpb) and SW (green) 
aides. The errors correspond to tha rariduab of the Gaua- 
rian fits divided by the ‘throughput” far the dbk FWHM 
as determined in the Appendix On both sides, the dbk 
FWHM steadily decreases from ~ 0.25" at r - 1.5" to 
• 0.15" at r - 1". Interior to r - l", the dbk FWHM 
fluctuates about a constant value, though this behavior 
m likely due to the difficulty of fitt ing a Gaussian pro- 
file in speckle-dominated regions, not bona fide clumpy 
structure. Beyond r = 1.5", where the dbk emission ap- 
proaches the photon noise limit, the FWHM estimates 
fluctuate wildly. 

3.2. Duk Position Angle 

lb quantify the “bow” structure easily seen in Fig 
ure 1 and identify any additional warping, we calculate 
the disk position angle aa a function of stellocentric d» 
taace for both the NE and SW sides using the dbk mid- 
plane pixel locations from the Gaussian fits described in 
§3 1. The position angle uncertainty al each radius re- 
sults from the difference between tbs Gaimsian fitted dbk 
midplane location at the radius and a m>d plane location 
defined by the brightest pixeL Here we formally assume a 
systematic uncertainty of 0.009* from the Yelda et al. as- 
trometric calibration, although fitting errors always dom- 
inate 

Beyond r = 0.9" , both skies of the dbk maintain a con- 
stant position angle, although they arc misaligned by 3-4 
degrees. Between 0.3" and 0.9", though, the disk emis- 
sion on both sides curves toward* the north, with the 
position angle decreasing on tho NE aide by more than 
20° and increasing by «v 10* on the SW aide. On the 
SW aide, this curvature b not continuous, leveling off at 
232* at 0.5-0 8" before resuming al smaller separations 
Beyond r = 1.6", the disk exhibits no obvious curva- 
ture/warping, though photon noise degrades the preci- 
sion of our estimates 

31 Disk Surface Brightness Profile 

lb calculate the disk surface brightness (SB), we fol- 
low Rodigaa et al (2012) and determine the median sur- 
face brightness in mJy foresee 3 in a 18-pixel diameter 
circular aperture We determine uocmaiutiee in the me- 
dian dbk surface brightness (also in mJyfarcoex *) at 
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a g.vrn angular separation la a way analogous to that 
wh *Ji are used to determine the disk SNR /pixel Specif 
ical.y, at each ptxel radius corresponding to a Gaussian 
fitted disk mid plane location we calculate the uncertainty 
in SB within non-overlapping circular apertures ldenti 
cal in else to that which we use to determine the disk SB 
and ooverlng all azimuthal angles exterior to the disk. 
We adopt the standard deviation of these SB estimates 
sa the uncertainty in the disk SB at each radius 
Figure 4 displays the surface brightness profiles for 
bota sides of the disk On both sides, the disk steadily 
brightens from r » 2.5" to r =■ 1". However, at smaller 
separations the profiles change dramatically On the NE 
side the disk has a roughly constant brightness at 0.5- 
0.9' before brightening from 3 to 6 mJu/oreeer 3 at 0 5" 
to G.4". The SW side of the disk displays a similar be- 
havior: a nearly constant brightness at 0 7-0 9", a possi- 
ble alight dip in brightness at 0.7", sad a sharp jump m 
brightness by a (actor of 4 &om 0.7* to 0 4*. Thus, tbs 
disk SB profile appears ‘wavy* intenor torsi*, 

The two sides of the dick differ slightly m toms other 
important respects The locations of the poaks/piateaus 
closer to the star for the NE side than for the SW side 
Tbs shifted profiles are consistent with the dust ring be- 
ing located at different stellocemnc distances between 
the NE and SW sides: a pencenter offset. The NE side 
also plateaus and may drop slightly at r ~ 0.3-0.35". 
Conversely, on the SW side the disk continues to brighten 
nil tne way to our inner working angle of 0.3" 

Moreover, our analysis confirms evidence for a bright- 
ness asymmetry between the NE and SW sides. These 
differences are significant even at small separations where 
we cetect the diak at a lower SNR. Interior to 0.5", the 
mad .an uncertainty m the disk SB is ~ 2.8 on the NE 
•ids and 4.6 on the SW side While these uncertainties 
are large, they are significantly smaller than the bright 
ness differences «tf> 0.5-0 6* where the SW’ aids is 2 3 
times brighter. At these angular sep a r ations, the ■*■/- \-c 
range in SB between tbs two sides do not overlap The 
SW side of the diak Is 50% brighter at r - 0.3-0. 4*. al- 
though the large error bars for the NE fade SB make this 
difference leas statistically significant 
We model the surface brightness profiles over r » 0.3- 
2.5" to derive power law indices a and 6 assuming a func- 
tional form of f(X) = aXV As a first guess, we take 
the logarithm of the surface brightness, adopt uniform 
weighting to each datapoint, and estimate the power in 
dicai by fitting a straight line and deriving the slopo. Wo 
then refine our power law index estimates by performing 
a Lr-enberg-Marquardt minimization, while considering 
measurement errors. 7b more precisely estimate the lo- 
cations of the power law breaks, we perform the above 
steps iteratively, varying the break locations and adopt- 
ing ones that minim ae the reduced g* between the data 
and the model. 

Table 1 summarises our results We (ail to identify 
any power law to either lobe that quantitatively fits the 
entire radial extent of the disk Fbr the NE (SW) side, 
the goodnesa-of-fit value (R 7 ) declines to aero at r m — 
10" 0.9"). Assuming a single poswr law for both sides, 
we derive best fit coefficients of a * 2.13, b — *6.01 at r = 
1 .05-2.6" for the NE side and a ~ 1.72. b - -5.49 for the 
SW' side at r ■- 0.95-2 5". However, the goodness of- fit 


i nterior for the NE side is low (R* ~ 0 75), indicating 
that at least this side may be beat modeled as a broken 
power law (see also Schneider et al 2005; Debea et al 
2009). 

For the NE side, the best-fit indices assuming a broken 
power law are a * 2.23, b - -6.19 at r m 1 1-1.4" and 
o . 1.31, b = -5.13 at r - 1.4-2 5" Fbr the SW side, 
the best-fit indices are similar a « 1 78, 6 w -5.71 at 
r - 1-1.6" and o =» 1.66, b - -5.33 at r =■ l 65-2 5" u . 
Formally, our values for b have a large uncertainty ainr* 
the goodncs»-of-fit criterion remains above 0.95 for a ± 
20* and 5 ± 0.5. 

We can compare our SB profiles to those from 
Schneider et al. (2005) and Debea at al (2009) to under- 
stand how the profiles change with wavelength In gen- 
eral, the shape of the SB profiles show good ag r ee m e n t 
with those derived from 1.1 fsm data by Schneider et ah 
(2005). Although they do not draw attention to any 
steep increase in SB at r * 0 5- 0.7*, their Figure 2 pro- 
vides some evidence for this feature, at least on the SW 
aide While we do not find evidence for a steep drop in 
SB at r = 03-0 5" for the NE aide as they do, the SB 
does plateau at r m 0 4" and drop slightly. The 2 05 /on 
profile from Debes et ah (2009) does not extend to r < 
0.5-0 6", so we cannot know whether they too find evi- 
dence for a jump in SB at small separations Our profiles 
are significantly steeper than those from Schneider et al. 
(2005) and Debes et al (2009) at r > 1 * obtained at 
shorter wavelengths (Figure 5). 

4. DEBRIS DISK MODELING 
41. Scattered Light Modeling 

To understand the disk’s grain properties and morphol- 
ogy, we compare our image to synethetic resolved disk 
models with a range al scattering properties We mod- 
eled a number density distribution of dust for the diak in 
the following cylindrical form, which allows for a variety 
of disk morphnbgMs: 

N(r,i)-«xp^^^^exp^^5^,r<r Wl (1) 

= (^) • r -’ w (J) 
where r a can be interpreted as tbs location of a birth ring 
of plane tesi main generating a collisionai cascade of dust 
with a characteristic width a. We assume a Gaussian 
scale height to the disk at all radii, and we allow for a 
power-law decrease in dust density beyond some radius 
n>r«*k. thus allowing for combinations of ring- like and 
disk like structures Wc place dust particles at various 
distances from the star following the density distribution 
above and distribute them uniformly in the azimuthal 
angle (9). 7b derive the integrated surface brightness, 
we project the density distribution on a two dimensional 
plane after transforming the coordinates of dual particles 
to account far inclination and position angle and scale 
the result to the observed surface brightness of the disk 

While t be mart and outer Stud radii reported bsrs are t h» 
anm lor which a p aw i law 6t is nvoat applicable, we obtain similar 
results for slightly diflareal radii • (. br„ ■ 0 95" and r n , 

•in' 
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In the Keck images. The surface brightness, F, of each the measured profile appears to flatten at r — 0.4", espe- 

particle is determined from its distance from the star dally for the NE side. Furthermore, our scattered light 

and its scattering angle m assuming a two component model predicts that the disk emission should be roughly 

Henyey-G reenstdn phase function: axMymmetric However, at r < 0.6" the disk’s SW side is 

significantly brighter and the profile breaks appear offset. 

( , 2 The disk’s asymmetric SB profile breaks could indi- 

1 ] ' a * 9i + 1 ~ tfacate an Wmmetry in the disk’s distance from the star 

r 2 + x 2 ) [ *(1 + ff? - 2gi cost*/) 1 * 2 (1 + gf - dabft^fcftxoE of position angk: e g. a pericentcr offset 

Tb investigate whether a pericenter offset better repro 
A multi-component phase function for drcumstellar duces the SB profile, we reconstruct a scattered light 

dust may be favored. For example, the phase function model with identical grain scattering properties as our 

of zodiacal dust in the Solar System has been modeled two-component model but make the SW side of the disk 


with multiple Henyey-Greenstem components, includ- 
ing a significant hackscattcring component (Hong et al 
19851. Recently, observations of HR 4796A’s disk sur- 
face brightness as a function of scattering angle showed 
an incredibly fiat inferred phase function for its dust at 
scattering angles > 50° (Dalle Ore et al. 2011), and ob- 
servations of the protoplanetary disk HD 100546 at mul- 
tiple wavelengths require forward scattering grains where 
the phase function becomes flatter at scattering angles 
> 38° (Mulders et al. 2012). 

Figure 6 displays the model that best reproduces the 
overall disk morphology and two rejected models. We 
attempted several possible structures and combinations 
of parameters, obtaining a good fit to the data with the 
model presented in Table 2 and shown in the left panel. 
A more rigorous exploration of parameter space and their 
possible degeneracies is beyond the scope of this paper 

Tb account for its warped, bow shaped appearance, the 
disk must contain some highly forward scattering grains 
at r ~ 110 AU, which cause a brightness asymmetry be- 
tween the front and back sides of the disk. At low scat- 
tering angles (small projected separations) the brightness 
asymmetry Is more pronounced than at larger projected 
separations, causing a change in the midplanc position 
angle. This type of “warping” structure is also seen 
in HD 15115 (Debesetal 2008; Rodigaa et al. 2012). 
However, the sharp breaks in the surface brightness at 
~1" are hard to reproduce with typical single component 
Henyey-Grcenstein phase functions (middle panel), re- 
quiring a flatter phase function at larger scattering angles 
(left panel). The surface brightness breaks could con- 
ceivably be reproduced by a superposition of two rings 
of isotropically scattering dust (right panel), but such 
a configuration does not give rise to the disk's warped 
appearance We cannot exclude the presence of a sec- 
ond inner ring (l.e. at r ~ 35-50 AU) for our modeled 
dust phase function, especially if the dust in the diak at 
wider separations is slightly less forward scattering than 
we have modeled This would allow the inner disk to 
dominate the surface brightness at small projected sep* 
rations. 

Although our two-component forward scattering model 
Is simple, it reproduces some key disk features. Figure 7 
and 8 (left panel) compares this model to the observed 
SB profile. The model clearly succeeds in reproducing 
the disk warping at r < 1" and the surface brightness 
profile break/turnover at ~ 110 AU. The model grains’ 
strong forward scattering at small angles causes the disk 
to appear very bright again at 6mall projected separa- 
tions (r < 0.6"). This feature agrees with our measured 
SB profile (Figure 8), though taken at face value the 
model predicts that this brightness accelerates, whereas 


5-10 AU closer to the star than the NE side (Figure 8). 
This model predicts that the NE side’s SB profile break 
starts at wider separations (r ~ 1.1" vs. 0.9" for the 
SW side) and that the SW side should be brighter at 
0.5 0.9", in agreement with the observed SB profiles. 

The pericentcr offset model’s fidelity isn't perfect: the 
SW side is significantly underi ominous and the exact 
locations of the breaks ore not well reproduced. How- 
ever, its success in reproducing the asymmetric SB pro- 
file breaks motivates more detailed scattered light disk 
modeling In particular, all single ring models (regard- 
less of any pericentcr offset) predict that the disk’s sides 
should be equally luminous at r ~ 0 5-0.4". Our data 
indicate otherwise (the SW side is still brighter by ~ 
50%), suggesting the need for further modifications for 
our scattered light model. New, higher SNR images of 
HD 32297’s disk will clarify the angular separation range 
over which the SW side is brighter and thus provide im- 
portant constraints for future disk modeling. 

4.2. Spectral Energy Distribution Modeling 

For a separate but related probe of the HD 32297 disk 
properties, wo modeled the disk spectral energy distri- 
bution (SED) from point-source photometry including 
newly-available data from the Sjntzer Space Telescope, 
WISE mission, and AKAJU satellite (Werner et aL 
2004; Wright et al. 2010; Murakami et al. 2007) Table 
5 lists the photometric data we consider 

Our methods follow those outlined in Plavchan et a). 
(2009), where we identify the best-fit dust temper a- 
ture(8), grain properties, and effective emitting areas in- 
corporating a downhill simplex algorithm (“Amoeba”) as 
described in Press et al. (1992). We consider cases where 
the dust radia tes like a blackbody such that r 4HtUAU = 

(280K/T) a x y/L,u,/Lo and where the dust's emiasiv 
ity scales with the effective grain size beyond a criti- 
cal wavelength (Badunan et aJ. 1992), c ot A - ^, which 
can place grains of a given temperature at larger dis- 
tances. We solve for the disk model parameters that 
minimiz e the fit resid uals relative to the flux (rmsraj = 
y / (EILo Ai/Fluxi^/N). In all cases, we assume that 
the grain populations have a characteristic size: for sim- 
plicity, we do not consider a grain trize distribution. 

Even though we include photometry not previously 
modeled, our fits are likely to be highly degenerate (e.g. 
Manesa et al. 2008; Fitzgerald et al. 2007). Therefore, 
instead of Identifying the single best-fit model and 95% 
confidence interval, we use several separate SED model 
runs to explore more focused questions about the HD 
32297 disk properties: 

• Model 1 - A Single Dust Population - Here we 
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assume that only one dust population contr.butes 
to the disk* IR-to-mm emission and allow fi to 
vary. Our goal with this lit ia to determine whether 
HD 32297 'b disk must have more than one dust 
population. 


• Model 2 - Two Duat Populations, One Dust 
Bolt at a Fixed Location Here we include two 
dust populations and allow their emisehity law* to 
vary but fix them to the same location We place 
the beh at r » 85 AU, or roughly interior to where 
the disk surface brightness profile begins to flatten 
Our goal here ia to determine whether more than 
one dust location is required. 


• Models 3-5 - Two Dust Populations /Bolts, 
Fixed Kmissivity Hero we Identify the best- 
fitting model with two du*t populations in two dis- 
tinct bells, assuming either that the grains behave 
as black bodies or have a A -1 etniasivity law. This 
approach follows that of Maness et al (2008) and 
will allow us to assess whether their formalism pro 
vides a better match to the photometry than as- 
suming one dust populalioa/belt. 

• Model 6 - Two Dust Populations /Bella, 
Variable Emiasivity, Fixed Location Here, 
we fix the locations of tbs dust to the peak of the 
surface brightness profiles — 0.4" and ~ l" or 45 
AU and 1 10 AU, which roughly cover the locations 
of the SB plateaus in the NE and SW Bides Our 
goal here is to see whether we can identify a good- 
fitting disk model that identifies the local maxima 
in the surface brightness profile as the locations of 
two separate debris belts 


• Models 7-0 - Throe Dust Population* /Belt*. 
Variable Emiasivity, Variable/ Fixed Loca- 
tion - Maness et al (2008) suggest that three sep- 
arate duat populations are needed to fit the HD 
32297 SED. In one case, we fix all the belts to be 
at 85 AU and incorporate a third, wanner dust 
population varying in gralu site and emiauvity. In 
two cases, we fix the outer two dust belts at 45 AU 
and 110 AU, while varying grain size and emiasiv 
ity and incorporate a third warmer dust population 
varying in in radius, grain size, and emiasivity. 

Figure 9 display* some of our SED modeling results, 
which are reported in Table 4 For all best-fit models, 
the italla: effective temperature is T. ~ 7870 7890 K 
and the star has little extinction (Av < 0 03) Regard 
lens of our assumed particle emusmty law, a tingle dust 
population model (Model 1, top- left panel) poorly repro- 
duces the SED, *t was found previously (Fitzgerald et al 
2007), significantly underpredicting the flux density at 8 
22 pm. Formally, having two or three dust populations 
in one belt at 85 AU (l.e. top-right panel) significantly 
improves the fit (rmj r «i « 0.05-0 08). However, the re- 
quired grain sues are too small (*- 4 nm-0 2 pm) to be 
realistic. Thus, the dust likely arisen from more than one 
single- temperature dust population with given grain size 
and Kninsrrtty power law. 

Flu assuming blackbody grains but incorporating two 
dust components at different locations significantly im- 
prove the fit (rnu^ ** 0.10; top-right) compared to 


a single duat component model The fits imply that 
the dust responsible for 10 1000 ftm emission is at 1.2- 
22 AU or r -v 0.01-0.2", but Fitzgerald et al (2007), 
Moerchen et al. (2007), and Maness et al. (2008) show 
that tbe disk emission instead originates on scales more 
comparable to =s 30-100 AU. The broad range of or- 
bital radii is indicative of the significant model degen- 
eracies involved Assuming that al least one of tbe dust 
components has an emiseivity of 0 * 1 yields fits with 
larger grains (0.3-3 pm) but worsens the fit (rmi^ ** 
017 0-21) Thus, if tbs dust disk geometries approxi- 
mate thin, isothermal rings and c on sist e n t of grains dom- 
inated by a single, characteristic stae, the grains likely 
have emiasivity powrr laws somewhere between 0 and 1. 

The bottom panels of Figure 9 show that it might be 
possible for the debris rings to reside at locations equal to 
the scattered light brightness peaks (r m 45 AU, 110 AU) 
and reproduce the IR to mm disk emission. A two duat 
population model (Model 0; bottom-left) reproduces the 
SED well at 1-8 pm and 17 -1300 pm but underpredicta 
the disk emission at ~ 12 pm by — 50%. Starting with 
tbe Model 6 results and inc orpo rating a third dust pop 
ulaDon, we reproduce tbe SED if the dust is at — 14 
AU (171 K), although tbe «nws»vity power law for this 
population is unphyncal [fi ~ 7 3) 13 and the grain sizes 
for the coldest dust component are too small (Model 8). 
Howev er , allowing the emiasivity law and grain sizes for 
all three components to vary yields an excellent fit to 
the data (RMS ~ 0.00). The emiasivity laws and gram 
sixes for this model (Model 9; bottom-right) are reason- 
able [fi = 0.37-0.77; o *= 0.20 0.70 pm). Here, the third 
component is at — 1 AU with a temperature of as 430 
K and is responsible for tbs weak 8 pm excess and sub- 
stantially contributes to the disk emission at 12-17 pm. 
Thus, although it is not nectMory, it is at least poutkU 
to identify the surface brightness profile peaks at 45 AU 
and 1 10 AU with separate thin dust nags responsible for 
broadband emission at 8-1300 pm provided that there 
exists an unresolved, wanner dust component located in- 
terior to 45 AU. 

We emphasize that the model fits we report in Table 
4 are but some of many possihlle fits, the range of which 
is further limited by our input assumptions about the 
disk For example, because our model is set up only 
to consider infinitely thin, isothermal nogs, it is unclear 
whether the warmer dust emi s si on we identify originates 
in a separate ring or rather an annulus, tbe outer edge 
of which we see at under separations Given the ex- 
treme number of model degeneracies, the most we can 
my Is that 1) the debns emission must originate from 
more than one population at 2) multiple locations and 3) 
among the many possible configurations, the du«t pop- 
ulations may coincide with the surface brightness peaks, 
lb make further progress, we need high spatial resolu- 
tion, high SNR imaging of the disk at a wider range of 
wavelengths to derive much more stringent constraints on 
the dust siaea/scattering properties and location(s) of the 
durt population ,'s) (s.g. Debes et al. 2009; Rodigas et si 
2012 ). 

» Th» modal ran Swto J — 7 J m iha -ban-fit* modal bacatmr 
(te Amoeba coda tmu 3 aa a traa traa paruoetar. raganflam of 
whathai tha vshia is physical, and aspkXla a gap Us SED masrtza 

al — 10 to adjfM a bailor it 



6. COMPARISONS TO PREVIOUS HD S2297 RESULTS 

From optical coronagraphic imaging. K&lu et al 
(2006) find evidence for warping in the HD 32297’* disk, 
especially for the NR tide, where the dak at r > 900-000 
AU appear* swept back compared to mid plane regions 
at » mailer angular separations They attribute the warp 
to ISM sculpting of the disk ee it moves south through 
ISM material The warped structure we resolve curves 
in the opposite direction from that expected due to ISM 
sculpting (see images for HD 32297/15115/81006 from 
Debes et al 2009) and does so for both tides of the disk. 
F\u’ hermore, ISM sculpting should be prominent only at 
larger separations in our images, where the disk is weakly 
detected (r > 1 .5*), because the grains need thus to be- 
come entrained in the ISM flow (e g Debes et al. 2009). 
However, optical and near- 1 R imaging probe two very dif- 
ferent characteristic grain population*, and the else scales 
over which we restive the dak do not overlap Thus, our 
result* are not in conflict Rather, ISM eculpting may ex 
plait the optical Image but does not explain the near TR 
image. 

VVe confirm and clarify three major claims from previ- 
ous near-IR imaging (Schneider et aL 2005; Mawet at al 
2009, Debra et al. 2009). First, we verify the brightness 
aeymmetry found by Schneider et aL and Mawet et al. 
(2002) at small (r ~ 46-55 AU) separations and show 
that it persists down to r ~ 35 AU, albeit at a 
lower statistical significance. Second, we recover a 
break/deviation in the SB power law at r ~ 100 AU found 
by Debos et al. (2009). Our new results show that the 
SB profile at smaller separations does not follow an ex- 
tension of the power law describing the dkk at larger sep- 
arai.^ns instead, H appears wavy. Finally, Mawet ct al 
(2002) claim that the ME side of the disk emission is 
truncated at r ~ 0.6". We find that the NE SB profile 
flattens (not brightens) from r - 1* to r 0.5* As- 
suming that their coronagraph attenuates some flux at 
separations slighter greater than their inner working an- 
gle (0.4*), most of the NE side emission srtD be hidden 
from their view. In contrast, the SW side dearly begin* 
to brighten at wider separations (r ~ 0 6-0 6*) visible 
by Mawet et al. (2009). Thus, our results agree 

Thermal IR imaging from Fitzgerald et u (2007) and 
Mocrcben el al (2007) shows that the HD 32297 dkk 
has an Inner clearing devoid of grama slightly larger than 
those wo can probe with our data Fitzgerald et al. find 
that the disk exhibits a bilobed structure and most of 
its 11.2 pm ftniSBion originates from r ~ 0. 6-0.6*. char 
actor i* tic of a dust ring with an inner edge at w 60 AU. 
Mourchcn et al resolve the disk at 12 pm and 18 pm out 
to slightly wider separations (r < 13*). Based on the 
disk's brightness temperature, they likewise find evidence 
for an inner clearing, albeit one that is slightly larger (r ~ 
80 AU) The plateaus in our SB profiles could identify the 
inner boundaries of dust populations (e g. Rodiga* et si. 
2012) Under this interpretation, our image k consistent 
with a dust belt truncated at r - 1 * (» HO AU) and 
a second ooe at w 0. 4-0.6* (46-70 AU) Because these 
authors did not analyze their data is the same way it 
k difficult to compare their results between each other 
and their combined results against ours Still, their re- 
sults and ours are qualitatively consistent with the HD 
32297 dkk having at least one dust population truncated 
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al small separations 
In agreement with Maness et al (2008)’* mm study, 
we find evidence far at least two dust grain popula- 
tions, responsible for the mid IR excess etmation and fer- 
IR/suhmm excess, respectively (see also Mocrcben et al. 
2007) They are able to fit the SED from 1 to 100 pm 
but not 1300 pm, a discrepancy they explain by adding 
a third, coid dust population. Our modeling generally 
fits the SED from 17 pm to 1300 pm but slightly under- 
predicts emission at 12 pm, a discrepancy we can solve if 
there exists an additional, unresolved warm dust popu- 
lation Our different results are likely a byproduct of our 
modeling assumptions they allow the dust populations 
to be spatially extended but fix the partick cmkaivity 
Iaw, whereas we assume the dust k confined in kotber 
mal rings but allow the emjsnvtty law to vary 
As noted in $3, the peak of HD 32297’s tnm emission 
from Manes* et al lies close to ths bright inner disk re- 
gion on the SW sufo. While formally the position *n g> 
(**» 46*) of the mm emission k offset by * 10 degrees, 
it tracks the disk's position angle at wider separations 
much better The disk’s warped appearance st r < 1 * k 
due to its strong forward scattering at small angles. If 
the disk’s grains instead isotropically scattered starlight, 
it# K. emission would lie almost directly on top of the 
mm peak Thus, we identify the mm and near IR bright 
ness asymmetry as originating from the same location 
Assuming the grains responsible for both the mm and 
near IR emission are likely the result of collisions, they 
may trace the same parent population of planetesimak. 

5.1. Summary and Future Work 

L sing Keck/NTRC2 K, coronagraphic imaging, we re- 
solve the HD 32297 dkk at a high «gnoMo- noise from 
r * 0 y' to r * 2.5*. We determine basic prop- 
erties (SH, FWHM, position angle), compare our image 
to dkk models with a range of (grain) scattering proper- 
tie#, and model newly -available, broadband photometry 
to provide a complementary investigation of HD 32297 e 
circumstellar environment. We obtain the following ma- 
jor results: 

• 1. We discover that HD 32297’s debris dkk ex- 
hibits a prominent warped or “bow" -shaped struc- 
ture interior to r ~ 1 * (~ 110 AU). 

• 2. Our new data clarifies the disk’s surface bright- 
ness profile at small separations Wo find that it 
has a “wavy" profile interior to r 1 10 AU with & 
plateau extending to r ~ 0. 5-0.7" (66-80 AU) be- 
fore the dkk brightens by factors of 3-4 at smaller 
separations. 

• 3. The dkk exhibits significant asymmetries be- 
tween the two side* The SW tide k brighter at 
r ~ 03-0. 6* by 50-150%, with the most statisti- 
cal differences being at r ~ 0. 5-0.6* Then* sepa 
rations are roughly where ths disk s mm emuaion 
peaks and are consistent with previous results (e g 
S chneid e r et al 2005). Our new analysis tiWifi — 
new asymmetnea. revealing that the locations of 
the peaks/ plateaus in surface brightness are likely 
shifted between the two sides, consistent with non- 
aximuthally symmetric structure. 
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• 4. A disk model with a fiat phase function and 
strongly forward scattering grains where the dust 
ring is centered on the star reproduces the "bow" 
structure, marginally reproduces the “wavy" SB 
profile and fails to reproduce the NE/SW asym- 
metries. A disk model with a 6-10 AU pen center 
oflset reproduces the asymmetric SB profile breaks 
interior to r < l", although its brightness asym- 
metry to limned to 0.5-0 9" and its match to the 
SB at wider separations is far poorer Thus, dust 
scattering plays a critical role in explaining key ob- 
served disk properties, but it l* unclear whether it 
explains all of the disk's properties we identify. 

• 5. HD 32297 must be surrounded by more than 
one dust population likely arts mg from different lo- 
cations in the disk. Although these populations 
need not be identifiable from our image, we can fit 
the diak SED by placing dust populations at the 
locations of the SB peaks from our image provided 
that there exists additional warm dust that we can- 
not yet resolve 

• 6 . The diak's brightness peak at r a* 0.4" coincides 
with the peak mxr. emission (Maaeas at al. 3008) 
If the grains responsible for both peaks are the re- 
sult of collisions, they may trace the same parent 
population of planeteumals 

It. summary, the HD 32297 disk appearance is 
broadly shaped by ISM interactions (Kalas at al 2005; 
Debes it al. 2009) at wide separations and by its grain 
scattering properties at small separations (this work). 
However. K is unclear whether either of these features 
by themselves explain the disk’s SB profile and (espe- 
cially) the disk's brightness asymmetry Furthermore, 
SED modeling provides evuience for mulUp> dust popu- 
lations, possibly multiple behs. lb explain the diak's SB 
profiie/tnultiple dust populations and (especially) bright- 
ness asymmetry, we may need to invoke other mecha- 
nisms 

Planets can sculpt dust into debris rings (# 4 . Kales 
200S. Quillen et al. 2006; Kalas et al. 2008) Further- 
more, planets can trap dust into resonant structures, 
which can appear as bright, overdenae regions like the 
brightness peak imaged here and in the mm at r — 0.4" 
(e.g IJou and Zook 1999; Kuchner and Holman 2003; 
Stax/, and Kuchner 2008). As argued by Wyatt et al. 
(2006) and Man ess et al (2008), the detectability of the 
resonant structure may be wavelength dependent. Large 
grains producing mm emission are fragments of the col 
tiding planeteaimals in resonance, but cannot be rapidly 
removed by radiative forces Thus they can trace the par- 
ent body resonant structure. Small, (sub)- micron sued 
grains likewise may trace resonant structure since they 
are preferentially produced in the most ooUtotonaily sc 
live, highest density regions and are otherwise quickly 
removed by radiation pressure. Grains with ustermedi- 
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ale sins producing emission at intermediate wavelengths 
(e.g. thermal IR) can be pushed out of resonance by ra- 
diation pressure /PR drag but are not stnall enough to 
be rapidly removed from the system (Wyatt et al. 2006, 
though see Kuchner and Stark 2010) Thus, it is possible 
that our near-IR image, when combined with the mm 
image, dentifiea planet- induced structure. 

Further near-IR and (sub)mm imaging is required 
to verify whether the SB profiles and brightness 
asymmetries are bona fide evidence for an embedded 
planet The current state-of-the-art nearto-mid IR 
high-contrast imaging facilities like the Large Bmocu 
lar Telescope have already shown great promise for re- 
solving disk's like HD 3 2297' • in scattered light and 
imaging self-luminous gas giant planets (Rodigas et al 
2012, Skemer et al. 2012) and upcoming planet im- 
agers l Jut SCEzAO on Subaru, GPI on Gsnum- South, 
and SPHERE on the VLT will be even more capa- 
ble (Esposito et al 2011; Martioache and Guyon 2009; 
Macintosh et al. 2008, Beuxit et al. 2008). Moreover, the 
thermal IR is well suited for imaging cxoplanets with a 
wide range of ages (e.g. Maroto et al. 2010; Currie et al. 
201 lb; Rodyas et al. 201 1), especially for stars like HD 
32297 whose bright, edge-on debris d.ak degrades planet 
sensitivity limits in the near IR Imaging with the Ata 
coma Large Millimeter Array (ALMA) can resolve the 
HD 32297 debris disk up to (actor of ~ 50 better than 
the CARMA observations reported by Maneas et al. and 
thus will provide a for better probe of any planet induced 
structure in the HD 32297 debris disk 
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TABLE I 

SiWAcm BwauTNEsa Profile Powkr Law 


iNoicni 

Angular Separation (") 

a,b(NE) 

a.b(SW) 

Single Power Law 
r - 1-2 5 

2.13, -6.01 

1.72, -5.49 

Broken Power Law 
r - 1.1-14 

r = 1-1,8 

2 23,6 19 

1.78, -5.71 

r = 1.4-2 5 
r = 1 65-25 

131. -6. 13 

1.66,-5 33 


Note. — Surface brightr.esi. profiles foe the 
NE sod SW ttdm determined by * Levcnbrrg. 
Marquardt leaat-squarea fit to the data. We assume 
a simple power law functional form of J{X) - aX*. 
Powr laws ta agrenment with tbs data at the 95% 
confidence limit include values for a that vary by ~ 
10% and values tor 6 that vary by ~ 0 5 dax. 


TABLE 2 


HD 33287 Scattered Light Model Parameters 


F# 

no AU 

«% 

13 AU 

KWHM, 

30 AU 


2 AU 

FWHM. 

5 AU 

Fkeeak 

110 AU 

r. 

096 

A; 

09 

n 

-0.1 

oa 

0.1 


3.4 

« 

r 

Poaition Angle 

47 5° - 


Currie. T., FVihagawa, M., at al., 2012, ApJ submitted, 
or XI v: 1206.0483 

Dalle Ore, CM. et al . 3011, AAA, 533, 96 
Dawson, R,, Murray-Clay, R_, Fabry cky, D., 2011, ApJ, 743, 17L 
Debw. J H , et al., 2006, ApJ, 664. 41L 
Debss. J. H , at al.. 3009, ApJ. 703. 318 
Esposito, S , et al.. 2011, SPIE, 8149, 1 
Fits nr aid, M . st al., 2007, ApJ, 670. 557 
Fitagerald, M., et aL, 2010, Proneedings of the conference In the 
Spirit of Lyot tOlO : Direct Detection of EaoploneU and 
CvTurruletlar Diets, October 26 - 29, 2010, University of Paris 
Diderot, Paris, Franc* Ed Anthony RoccalettL 
GoUmowski, D., et al., 2006, AJ, 131, 3100 
Heap, S.. et al.. 2000, ApJ, 539, 436 
Hines, D., et al.. 2007, ApJ, 671, 1661. 

Hong, 8. 8., 1965, AAA. 146, 67 

Kalaa, P., Graham, J., and Clam pin, M., 2006, Nature, 435, IU67 
Kalas, P. 2005b, ApJ, 635. 1601. 

Kalaa, P., et al 2007, ApJ, MI, SSL 
Kalaa, P , at aL, 2008, Science, 322, 1345 
Kenyon, 3. J.. Bromley. B , 2008, ApJS. 170, 451 
Kocfcnar. M., Holman. M , 2003, ApJ, 568, 1U0 
Kudaer, M., Stark, C . 3010, AJ, 140, 1007 
Lafrctuare, D., et al., 2007, ApJ, 660, 770 
Lupine*, AM. et al , 2010, Sdence, 329, 57 
Llou, J. C., Zook, H., 1999, AJ, 118, 580 
Macintosh, 3 . et al . 2008, SPIE, 7015, 31 
Martinadte, F . Guyon, O , 2009, SPIE, 7440. 30 
Mantas, H , et aL, 2008. ApJ, 686, 25L 


Marc**, C , et aL, 2006. ApJ, 641. 556 
Marota, C . st al., 2006, Science, 322, 1348 
Marob, C., et al., 3010, Nature, 466. 1080 
Mawet, D, et al., 2009. 702, 47L 

Metchev, S„ Maras. C . Hillenbrand, L., 2009, ApJS, 181, 62 

Moerchen, M., et al , 2007, ApJ, 666, 109L 

Mulders, G-. et al., 2012, AAA submitted 

Murakami. H., et al.. 2007, PASJ, 50, 369 

P'.avchan, P., et a L. 3009, ApJS, 608, 1068 

Presa, W., et al., 1902, Numerical Retypes m Fortran : The Art of 
Scientific Computing. Cambridge: University Preas, 1992, 2nd 
ed 

Quiflen. A , 3006, MNRAS, 372, 14L 
Rhee, J., et aL, 2007, ApJ, 660. 1556 
Rodigaa, T. J.. et al.. 2011, ApJ, 732, 32L 
Rodlgaa, T J.. et aL. 2012, ApJ la proa 
Schneider, G., et al , 2005, 629, 117L 
Sanm, A., et al., 2012, ApJ submitted 
Smith, B-, Terrile, R- J., 1984. Sdence, 220, 1421 
Stark, C , Kuchner, M., 2008. ApJ, 686, 637 
Su, K., et al . 2000, ApJ, 705, 314 
Ttalmann. C.. et a l. 2011, ApJ, 743, 6L 
van Leeuwen, F, 2007. AAA, 474, 663 
Werner, M . et al., 2004, ApJS, 154, 1 
Wright. E L-, et al , 3010, AJ, 140, 1868 
Wyatt, M . et al., 1999, ApJ, 527, 918 
Wyatt, M-, et a L, 2006, ApJ, 839. 1153 
Wyatt, M., 2000, ARAAA, 46. 330 
Yalda, S., et al.. 2010, ApJ, 725, 331 



TABLE 3 

HD 33397 Photometric Data 


Filter /Wavelength 

Flux (mJy) 

ffPlux (mJy) 

Source 

fl/0 4380 

1836 74 

28 75 

TYCHO-n'Uaaa) 

V/0.B45 

2012.91 

24 10 

TYCHO-Il(traaa) 

Jl\ 335 

134’. 80 

34 96 

2MASS 

///l 662 

91348 

50.58 

2MASS 

K.2159 

611.41 

11.95 

2 MASS 

3.37 

280.23 

8.21 

WISE 

4.62 

151.28 

3.78 

WISE 

8 

65 18 

275 

IRSep 

112 

49.9 

2.1 

Kltagerald <* al (2007) 

11 66 

53 

5.3 

Moeriiiec h al ( 2007) 

13 08 

55 25 

1 20 

WISE 

16 

71 28 

2 53 

IRSep 

18.3 

90 

13 6 

Moerrhen cl al (2007) 

22 19 

212 99 

5 55 

WISE 

24 

225.2 

482 

IRSep 

TO 

850 

60 

IRSep 

90 

823 2 

116 

AKAR1 

180 

<460 


IRSep 

1300 

5.1 

El 

OAHMA/Mra* « «I. (2008) 


Note. — TYCHO fj(tiwu) refer* to TYCHO- II catalog data transformed into 
the atandard Johnaon Cocaine photometric mtan IRSep refer* to the IRS en 
haoced product* dataaet ea queried from the NASA/IPAC Infrared S c ie nce Archive. 


TABLE 4 

HD 32397 SED Mooeunc Results 


Vfodel ID 

RMS 

fi-vt.: (AU) 



fl\ 

* 


aj (**m) 

<*3 

a* 

1 

0.23 

85 



0 94 

_ 

— 

0.005 



2 

008 

85 

85 

— 

0 72 

0.39 

— 

0004 

0.197 

- 

3 

0.16 

1 21 

21.66 


0 

0 

— 

039 

0.41 

— 

4 

a 17 

55 64 

527.56 


1 

1 

- 

0.30 

0.99 

— 

5 

0.21 

15 60 

97 50 

— 

0 

1 

— 

0.91 

286 

• 

6 

0 15 

45 

110 

— 

0.80 

0.31 

— 

0.39 

0 025 


7 

0.05 

85 

85 

85 

084 

0 61 

0.45 

0 006 

0024 

0.725 

8 

0.09 

14 

45 

110 

7.27 

0.80 

0 31 

1.5® 

0.39 

0.03 

9 

006 

1.1 

45 

110 

0.37 

077 

0.43 

0.70 

037 

0.20 


Not*. — RMS refer* to the fit residuals relative to the flux rma^, - ^(2jJLo A,/FJ u*,)*/W, 
where Af ia the number of flux density meaauremenu (18). R+,, t refer* to the duet ring'* itelloceatric 
dlatanoa, (i is the particle nmieaivity power law, and a a the grain alee In truer ans, where 3na « A* 
(»no Badunas at al 1993). Values in bold are fixed fur a given model run, whereae other* are ‘fitted’ 
value*. 
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Pic L— Raducad Image (lop) awl n|ntl-u> mar map (bottom) for otu NIRC2 HD 32207 data. The cokxbar depicting units for Urn 
imigi are la counta, wheraaa they range from 0 to 9o for the signal to- noise map The central dark region identifies the ooronagraphk: spot 
(r s 0.3"). Tbs panels have lha same alar scale 
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1‘1C. 2. (Top) Our Imago from the panel of Figure dLired tap lay«d withJIi different color dttetd > to belter Illustrate the r.gn-ftcanre 
of our dink detection at the smallest separations (r =0.3-0 6") and the brightnms asymmetry. (Bottom) Our image resampled to the 
•amt spatial resolution a a the Palasnar/WCS image from Mn«t et al. (2009). The green cress Ideal if.es the position (and positional 
uncertainties) of the peak brightness in the millimeter (Manem et al 2008) Both panels are displayed In units of counts 










Pro 0.- Scattered :.ght models incorporating diSerent grain properties/disk geometries: (left) a two- component Henyvy-Gremistnin 
modal with strongly forward scattering grains at small scattering angles bat weakly scattering ones at Larger angles (larger projected 
separations), (middle) a simple forward scattering grain model, (right) two isotropically scattering dust rings. The Lefthaod model best 
reproduces the disk SB profiles. The unit* are in counts 
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Fig 7. — Comparison* btlwta th« modeled end measured disk position angles (right). The triangle identifies the location ol the nur. 
bright saw peak from Manes* et aL (200S) We use the twxonponinl modal with the dust ring cantered on the star and with a 10 AO 
offset. 



FiG. Comparisons between the modeled and measured surface brightness profiles for a dust ring entered on the star (left) and ooe 
with a 5-10 AU per -center offset (right; the SW side is 5-10 AU doner). Both model* reproduce the wavy SB profile The offset makes the 
disk model qualitatively reproduce the differences In the observed SB peaks /plateau* at r < l", although tt degrades the model's fidelity 
at r > I ", especially an the SW aide, where it M substantially under turn ino us 






Fas 10 — Images used to model biasing tom LOCI processing. (Left) Our final imago with a fake disk added (Right) Our final image 
where we add the fake disk to each registered image and process the set of imager with our pipeline In both esses, we rotate the image to 
the parallactic angle of the first image in the sequence (PA ~ -15 08°), not to true north as ws do la P.gure 1. 














R (orc-*#c) 

Fio 12.— Difference In position angle for the (nice disk before end after processing Oar processing does not bine the disk astrometry by 
more than 1-2° at any separation. 
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APPENDIX 

CORRECTING FOR DISK PHOTOMETRIC/ ASTROMETRIC BIASES WITH LOCI 

The LOCI- based PSF subtraction approach can bias the photometry and astrometry of point sources like planets and 
extended structure like disks (Lafremere et al. 3007; Thalmann et al. 2011). For point source companions, imputing 
fake companions into registered images at a range of angular separations, processing tbeee images, and then comparing 
the output to input fluxes and centroid positions of the take sources corrects for these biases Correcting biases for 
disks not v iew e d perfectly edge-on /pencil- thin is more difficult since we do not know, a p ri o ri, the disk's true FWHM, 
our ignorance of which also limits our ability to ’ debts* ' other disk properties whose determination depends on the disk 
HWHM (Le surface brightness). Here we describe our method for mitigating LOCI- based photometric/ astrometric 
biases, which largely follows that of Rodigaa et al. (3012). 

We first construct model disks each with midpUae brightnesses (and thus signal- to nowe ratios) about twice that of 
the pro ces sed, real disk image and surface brightness profiles comparable to the real disk profiles beyond r =* 1*. We 
comirier two FWHM values of 0.2" and 0 25" at 1", or ~ 20 10% larger than the values we get between r = 0.3 and 
1* jrlor to applying the corrections we derive in this section (Figure 11, top panels). F\»r both runs, we imput the fake 
disk into registered images, rerun our processing pipeline and then compare the output and input surface brightness 
profiles, disk thicknesses, and disk position angles. 

Ai shown in Figure 1 1 (middle and bottom panels), our processing mlnimially biases the thinner model disk (FWHM 
■ 0 2") ex t erior to r = 1 , reducing its surface bright new and thickness by no more than 20% Interior to r ■ 1", the 
surface brightness and thickness drop to no less than GO-75% of their original values Biasing for the FWHM * 0.25" 
disk (not shown) is only slightly more severe, resulting in an additional ~ 5% (10%) drop al r > (r < 1") in SB 
and FWHM. Fbr both fake disks, LOCI minimally biases the disk position angle measurements (Figure 12). 

The ‘observed” FWHM for the model with the initially thinner disk (FWHM =0 TT) Is ~ 10% smaller than the 
‘observed FWHM for the initially thicker disk (FWHM * 0.25") and more similar to that which we find for the rnal 
HD 32297 disk prior to applying any bias correction Thus, we derive throughput correction* for the disk surface 
brightness and FWHM for the thinner disk, fitting a simple, unweighted power law to data betwe en r = 0.3" and r ■■ 

The deviation in throughput for individual points vs, our throughput function are u largo ss 10-20%. However, 
our modeling errors for the disk SB and FWHM are larger in regions where biasing hum LOCI is important and at 
r > 1.5-2" where the disk Is intrinsically much fainter Therefore, we leave s more detailed, robust calibration of 
disk parameters from LOCI -processed images to a future work where a higher SNR disk detection at r < 1" will help 
improve our fitting precision. 



